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Using a Bayesian method for solving inverse quantum problems, potentials of quantum systems 
are reconstructed from coordinate measurements in non-stationary states. The approach is based 
on two basic inputs: 1. a likelihood model, providing the probabilistic description of the measure- 
ment process as given by the axioms of quantum mechanics, and 2. additional a priori information 
implemented in form of stochastic processes over potentials. 
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The first step to be done when applying quantum me- 
chanics to a real world system is the reconstruction of its 
Hamiltonian from observational data. Such a reconstruc- 
tion, also known as inverse problem, constitutes a typical 
example of empirical learning. Whereas the determina- 
tion of potentials from spectral and from scattering data 
has been studied in much detail in inverse spectral and 
inverse scattering theory jlJB], this Paper describes the 
reconstruction of potentials by measuring particle posi- 
tions in coordinate space for finite quantum systems in 
time-dependent states. The presented method can easily 
be generalized to other forms of observational data. 

In the last years much effort has been devoted to many 
other practical empirical learning problems, including, 
just to name a few, prediction of financial time series, 
medical diagnosis, and image or speech recognition. This 
also lead to a variety of new learning algorithms, which 
should in principle also be applicable to inverse quan- 
tum problems. In particular, this Paper shows how the 
Bayesian framework || can be applied to solve problems 
of inverse time-dependent quantum mechanics (ITDQ). 
The presented method generalizes a recently introduced 
approach for stationary quantum systems QJ^]. Com- 
pared to stationary inverse problems, the observational 
data in time-dependent problems are more indirectly re- 
lated to potentials, making them in general more difficult 
to solve. 

Specifically, we will study the following type of obser- 
vational data: Preparing a particle in an eigenstate of the 
position operator with coordinates xq at time to, we let 
this state evolve in time according to the rules of quan- 
tum mechanics and measure its new position at time t\, 
finding a value x±. Continuing from this measured posi- 
tion xt, we measure the particle position again at time 
t-2, and repeat this procedure until n data points X4 at 
times ti have been collected. We thus end up with obser- 
vational data of the form D ={(xi, A,-,a;,'_i)|l < i < n}, 
where Xi is the result of the i-th coordinate measurement, 
A; = ti~ ti_i the time interval between two subsequent 
measurements and Xi-i the coordinates of the previous 
observation (or preparation) at time 

We will discuss in particular systems with time 
independent Hamiltonians of the form H = T + V, con- 



sisting of a standard kinetic energy term T and a local 
potential V(x,x') — S(x — x')v{x), with x denoting the 
position of the particle. In that case, the aim is the re- 
construction of the function v(x) from observational data 
D. (The restriction to local potentials simplifies the nu- 
merical calculations. Nonlocal Hamiltonians can be re- 
constructed similarly.) 

Setting up a Bayesian model requires the definition of 
two probabilities: 1. the probability p(D\v) to measure 
data D given potential v, which, for D considered fixed, 
is also known as the likelihood of v, and 2. a prior prob- 
ability p(v) implementing available o priori information 
concerning the potential to be reconstructed. 

Referring to a maximum a posteriori approximation 
(MAP) we understand those potentials v to be solutions 
of the reconstruction problem, which maximize p(v\D), 
i.e., the posterior probability of v given all available data 
D. The basic relation is then Bayes' theorem, according 
to which p(v\D) cx p(D\v)p(v) . 

One possibility is to choose a parametric ansatz for 
the potential v. In that case, an additional prior term 
p{v) is often not included (so the MAP becomes a max- 
imum likelihood approximation). In the following, we 
concentrate on nonparametric approaches, which are less 
restrictive compared to their parametric counterparts. 
Their large flexibility, however, makes it essential to in- 
clude (nonuniform) priors. Corresponding nonparamet- 
ric priors are formulated explicitly in terms of the func- 
tion v(x) Q. Indeed, nonparametric priors are well 
known from applications to regression j?], classification 
H , general density estimation |)| , and stationary inverse 
quantum problems 0,0] . It is the likelihood model, dis- 
cussed next, which is specific for ITDQ. 

According to the axioms of quantum mechanics the 
probability that a particle is found at position xi at time 
ti, provided the particle has been at at time tj-i, is 
given by 



Pi = p{x i \&. i ,x l -i 1 v) = \(j>i{xi)Y 
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where 



<j>i{xi) =< Xi \4>i> = <Xi\UiXi-x>, (2) 
are matrix elements of the time evolution operator 
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setting h — 1. The transition amplitudes can be cal- 
culated by inserting orthonormalizcd eigenstates ipa of 
H, with energies E a , 



4>i(xi) 



-iAiE c 



(4) 



Clearly, it is straightforward to modify (Q) for measur- 
ing observablcs different from the particle position. It 
is also interesting to note that the transition probabil- 
ities (0) define a Markoff process with Wi(x — > x') — 
p(x'\Ai, x, v). For real eigenfunctions ip a (x), i.e., for a 
real Hamiltonian with real boundary conditions, they 
obey the relation Wt(x — > x') = Wi(x' — > x). It follows 
that the detailed balance condition, p s tat{x)Wi(x — > x') 
= Psta,t{x')Wi(x' is fulfilled for a uniform p sta t(^), 

which therefore represents the stationary state of the 
Markoff process of repeated position measurements. 

Having defined the likelihood model of ITDQ, in the 
next step a prior for v has to be chosen. A convenient 
nonparametric prior p(v) is a Gaussian 
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with (real symmetric, positive semi-definite) inverse co- 
variance Ko, acting in the space of potentials, and mean 
vo(x), which can be considered as a reference potential 
for v. Typical examples are smoothness constraints on 
v which correspond to choosing differential operators for 
Ko- Reference potentials can be made more flexible by 
allowing parameterized families vq(x; 9). Within the con- 
text of Bayesian statistics such additional parameters 9 
are known as hyperparameters. In MAP approximation 
the optimal hyperparameters are determined by maxi- 
mizing the posterior (Q) simultaneously with respect to 
9 and v(x) ||To|| . A simplified procedure consists in us- 
ing a parametric approximation v(9) which maximizes 
the likelihood IliPj(^) as reference potential vo for the 
nonparametric reconstruction v(x) (51. 

If available, it is useful to include some information 
about the ground state energy Eo(v), which helps to de- 
termine the depth of the potential. This can, for exam- 
ple, be a noisy measurement of the ground state energy 
which, assuming Gaussian noise, is implemented by 



p E CX e~t (Ea(v)-.f 



(6) 



Combining (||) and ^ with ([!]) for n repeated coordi- 
nate measurements starting from an initial position xq, 
we obtain for the posterior ([?]), 



p(v\D) cx pg(v)pe(v) Y\_Pi 

i=l 



(7) 



To calculate the MAP solution v* —a,rgmax v p(v\D) we 
set the functional derivative of the posterior (Q), or tech- 
nically more convenient of its logarithm, with respect to 
V, denoted S v , to zero. This yields, 

= S v ]np(v\D) = S v \np G (v) + S v hip E (v) + ^ S v lnpj, 



(8) 



with 



6 v lnp G (v) = -K (v - wo), (9) 
S v lnp E (v) = -fi(E (v) - k) S v E (v), (10) 
5 v \npi = 2Re[(j)^ 1 {xi)5 v 4>i{x i )]. (11) 

The functional derivative 5 v 4>i can, according to Eq. (0), 
be obtained from S v ip a - The still required 8 v ^ a and S v E a 
can then be found by calculating the functional derivative 
of the eigenvalue equation Hip a — E a ip a . Using 



5 v{x) V(x', x") = 8{x - x')8(x' - x"), 



(12) 



5 V ( X \ denoting the x component of functional derivative 
S v , we find, 



8 v (x)E a = <ip a \ 5 V ( X )H \ip a >= \tp a {x)Y 
8v(x)ip a (x') = ^2 - ip*y(x')ip*(x)il> a (x). 



(13) 
(14) 



Collecting the results, gives 

8 v (x)4>i{xi) = 5 v (x) <x i \U i x i - 1 > = 

(-iAi\ip a (x)\ 2 )) ip a (xi)^(xi-i) 
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1p*(Xi-l)lp 7 (x)lp^(x)tp a (Xi) 



(15) 



Inserting Eq. (]13j) for a = in Eq. fll0|) and Eq. fll§) in 
Eq. ([ll]) a MAP solution for the potential v can be found 
by iterating the stationarity equation (Q) numerically on 
a lattice. Clearly, such a straightforward discretization 
can only be expected to work for a low-dimensional x 
variable. Higher dimensional systems usually require ad- 
ditional approximations pL 

As the next step, we want to check the numerical feasi- 
bility of a nonparametric reconstruction of the potential 
v for a one-dimensional quantum system. For that pur- 
pose, we choose a system with the true potential 



Vtrue(x) 



(16) 



/2trt 

where c\ = —10, c 2 = —2, and a = 2. An example of 
the time evolution of an unobserved particle in the po- 
tential Wtruo is shown in Fig. ^. As input for the recon- 
struction algorithm 50 data points X{ are sampled from 
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the corresponding true likelihoods p(a;|A,-, Xi-i, i>tr U e)- A 
corresponding path of an observed particle is shown in 
Fig.| 

Besides a noisy energy measurement of the form (^) we 
include a Gaussian prior (||) with a smoothness related 
inverse covariance 

3 2m / £)2 \ 

K ( X) x') = 5(x- X ')Xj2(-l) k ^- k ^ . (17) 

To obtain an adapted reference potential vq for the Gaus- 
sian prior, a parameterized potential of the form 



vo(a, b, c) — min[0, a(x — b) 2 + c], 



(18) 



is optimized with respect to a, 6, c by maximizing the 
"extended likelihood" J2i^ n Pi( v o) + ^Pe(vo)- Finally, 
the stationarity equation (|^) is solved by iterating ac- 
cording to 



v (r+i) _ v (r) _|_ 



w - w (r) + 



(19) 



K " 1 { 2 ErRe[<5 t , ln^(^)]+^(lnp G + lnp E )} 



The resulting nonparametric ITDQ solution witdq (see 
Fig. |^) , is a reasonable reconstruction of t>t rU c ) an d clearly 
better than the best parametric approximation vo(a, b, c). 
It is only the flat area near the right border where, due 
to missing and unrepresentative data, the reconstruction 
differs significantly from the true potential. 

Fig. ^ compares the sum over empirical transition 
probabilities — ~ x i) as derived from the ob- 

servational data D with the corresponding true ptruc = 
-tY^i=iP( x \Ai,Xi-i,v tr uc) and reconstructed pitdq = 
- Y^i=i P( x \A-i, Xi-i, ^itdq)- Due to the summation over 
data points with different Xi-±, the quantities shown in 
Fig. U do not present the complete information which 
is available to the algorithm. Hence, Fig. || depicts the 
corresponding quantities for a fixed Xi-\. In particular, 
Fig. compares the reconstructed transition probability 
(Q) with the corresponding empirical and true transition 
probabilities for a particle having been at time ti-% at 
position Xi-i —1. The ITDQ algorithm returns an ap- 
proximation for all such transition probabilities. 

Figs. [| and || show, that the reconstructed Uitdq 
tends to produce a better approximation of the empir- 
ical probabilities than the true potential wt rU c- Indeed, 
the error on the data or negative log-likelihood, £d(v) 
= — hipi(v), being a canonical error measure in den- 
sity estimation, is smaller for witdq than for ut rU o- A 
smaller A, i.e., a lower influence of the prior, produces 
a still smaller error £_d( w itdq)- At the same time, how- 
ever, the reconstructed potential becomes more wiggly 
for smaller A, being the symptom of the well known 
effect of "overfitting" . The (true) generalization error 
e g (v) — — Jdx dx'p(x)p(x'\x, Wtruc) liip(x'\x, v) [with uni- 
form p(x)], on the other hand, can never be smaller for 



the reconstructed witdq than for v tmc - As it is typical 
for most empirical learning problems, the generalization 
error £ 9 (witdq) shows a minimum as function of A. It is 
this minimum which gives the optimal value for A. Know- 
ledge of the true model allows in our case to calculate the 
generalization error exactly. If, as usual, the true model 
is not known, classical cross-validation Q and bootstrap 
JPl[ techniques can be used to approximate the general- 
ization error as function of A empirically. 

Alternatively to optimizing A or other hyperparame- 
ters one can integrate over them JTc| . Similarly, studying 
the feasibility of a Bayesian Monte Carlo approach, con- 
trasting the MAP approach of this paper, would certainly 
be interesting. 

In summary, this Paper has presented a method to 
solve inverse problems for time-dependent quantum sys- 
tems. The approach, based on a Bayesian framework, is 
able to handle quite general types of observational data. 
Numerical calculations proved to be feasible for a one 
dimensional model. 
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FIG. 1. Time evolution of an unobserved particle with mass 
m = 1 started at time to = from xo = in the potential (fl6|). 
Shown is the transition probability p(x\At = t, xq = 1, Vtrue)- 




FIG. 2. Time evolution of an observed particle with mass 
m = 1 in the potential (^). The figure shows for each data 
point i the probability p(x\Ai = 5, Xi-i, Wtrue), starting from 
xo = 0. (Hence, the probability at i —1 corresponds to that 
shown in Fig. hi at t = 5.) The actual data points X{ have 
been sampled from that probabilities and form the observed 
path shown on top as a thick line. 
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FIG. 3. Numerical reconstruction of a potential from 50 
coordinate measurements (see Fig. ^). Shown are the true 
potential utruc (thin line) , the best parametric approximation 
used as reference potential vo (dashed line), and the recon- 
structed potential «itdq (thick line). Parameters: m = 1, A; 
= 5, vtiuc of Eq. (|l6|), no of the form (|li|), Gaussian prior with 
Ko as in ( |l7| ) with A = 0.1 and cto = 3, pe (q) with /i = 10 
and k = -Eo(^true), periodic boundary conditions for i\) a , fixed 
boundary values v{— 10) = v(10) = 10 s for uitdq, vo and « trU c, 
calculated on a lattice with 21 points. Errors: £d(^itdq) = 
99.1, eo(wtrue) = 104.4, e g (v ITDQ ) = 1.891, e 9 (« truo ) = 1.818. 
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FIG. 4. Sum of empirical transition probabilities p e m P 
(bars), the corresponding true ptmc (thick line) for v tTUC , and 
the reconstructed pitdq (thin line). 
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FIG. 5. Same functions as in Fig. ^, but restricted to mea- 
surements of a particle which has been at position = 1 
at the time of the previous measurement. 
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